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Oceanic Kelvin and Rossby waves play an important role in tropical climate and El-Nino dynamics. 
Here we develop and apply a climate network approach to quantify the characteristics of El-Nino 
related oceanic waves, based on sea surface height satellite data. We associate the majority of 
dominant long distance (> 500 km) links of the network with (i) equatorial Kelvin waves, (ii) off- 
equatorial Rossby waves, and (iii) tropical instability waves. Notably, we find that the location of 
the hubs of out-going (~ 180°E) and in-coming 140°W) links of the climate network coincide 
with the locations of the Kelvin wave initiation and dissipation, respectively. We also find that 
this dissipation at ~ 140°W is much weaker during El-Nino times. Moreover, the hubs of the off- 
equatorial network coincide with the locations of westerly wind burst activity and high wind vorticity, 
two mechanisms that were associated with Rossby waves activity. The quantitative methodology 
and measures developed here can improve the understanding of El-Nino dynamics and possibly its 
prediction. 

PACS numbers: 


The El-Niho-Southern Oscillation (ENSO) is the 
largest climatic cycle on annual time scales, and is one of 
the most important processes that affect the natural cli¬ 
mate variability. El-Nino has a profound effect not only 
on the tropical atmospheric and oceanic conditions, but 
also on remote regions including North America, Indian 
monsoon region, and Antarctica [l}-{3] . Classical theories 
of El-Nino attribute its self-sustained dynamics mainly 
to positive and negative feedbacks between the depth of 
the thermocline, sea-surface temperature, and winds. 

Oceanic waves play an important role in El-Nino dy¬ 
namics, and the equatorial region can be regarded as a 
natural waveguide of several degree latitude width cen¬ 
tered at the equator ■ Q|. Tropical waves, such as Kelvin 
and Rossby waves [j], are propagating in this waveg¬ 
uide. The intraseasonal Kelvin waves have been exten¬ 
sively studied using various observational and modeling 
methods due of their importance to the theory of El- 
Nino @. For example, the 1997/8 strong El-Nino event 
was initiated by a group of very pronounced Kelvin waves 
which were generated by energetic westerly wind bursts 
(WWBs) Recently, satellite altimetry data has be¬ 
come an invaluable resource to study equatorial wave 
dynamics |T, 8j. The velocity of waves in the equato¬ 
rial Pacific region can be estimated using the so-called 
Hovmoller (longitude-time) diagrams, by transferring the 
data to the wavenumber-frequency space, or by decom¬ 
posing the data into several leading modes of equato¬ 
rial Kelvin and Rossby waves [&, . Generally speaking, 

the typical phase speeds of the first baroclinic mode of 
equatorial Kelvin and Rossby waves are ~ 2.8m/s and 
~ 0.9m/s, respectively Si where Kelvin waves exist on 
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the equator and propagate eastward while Rossby waves 
are off-equatorial waves that propagate to the west. 

Tropical oceanic waves connect different, sometimes 
very far, regions. The set of regions V and their in¬ 
teractions E can be encoded in a graph G{V , E) named 
in recent literature “climate network” 0 - Many fields 
of research benefit from a rich set of quantifications and 
algorithms that has been developed for analyzing such 
networks of interactions. Examples include social and 
biological systems, information flow through the world 
wide web and physiological activities (see [ill, 03 f° r 
reviews). In climate science, network theory has been 
used, e.g., to infer known climate patterns (like ocean 
currents), to investigate atmospheric variabilities, and 
to study low-frequency climate phenomena 0, I l/l 1 1 Of . 
It is common to regard links from E whose statistical 
significance is above a certain threshold as weighted or 
unweighted links. Since network based methods empha¬ 
size coordination between locations rather than local dy¬ 
namics, they exhibit an increased sensitivity to climate 
patterns which are not easily captured by conventional 
analysis [3, 01. In spite of the importance of oceanic 
dynamics on the climate system, only a few studies have 
used oceanic data within the framework climate net¬ 
works 0, 0| • Also, the studies of El-Nino using cli¬ 
mate network methodo logy took advantage of near sur¬ 
face atmospheric data 0, 0, 0] . Here we aim to fill 
this gap, focusing on the important ocean-atmosphere 
phenomenon, El-Nino. 

By transforming the daily and relatively long (~20 
years) satellite based altimeter (i.e., sea surface height) 
records into climate networks, we find that the time de¬ 
lays, the velocities and the directions associated with 
climate network links can be attributed to equatorial 
Kelvin, Rossby and tropical instability waves (TIW). As 
will be shown below, the topology of the climate network 
is strongly affected by El-Nino, which implies that the 
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relative number of links in the eastern equatorial Pacific 
is significantly higher during El-Niho times compared to 
normal times. Notably, the network properties reveal the 
locations of the wave initiation as well as its strong dis¬ 
sipation. Our results suggest that climate network may 
serve as a useful tool to study wave dynamics, and more 
importantly, network features may serve as their quanti¬ 
tative measurements. 

We analyze the daily sea surface height (SSH) anomaly 
data (2lj| . We focus on the tropical Pacific ocean be¬ 
tween 7°S to 7°N, and between 120°E to 280°E (80°W); 
the time extent of the data is from January 1 1993 to 
December 31 2010. The spatial resolution is 1° and 
4° in the meridional and zonal directions, respectively. 
Totally, there are 615 grid points, which are regarded 
as nodes of the climate network. For each node of the 
network, we analyze daily SSH values after subtracting 
the seasonal cycle. The filtered SSH record is denoted 
by H. We also define the normalized series Q s (t ) = 

[H s (t) - (H s (t))]/[(H s (t) - (H s (t))) 2 ) 1/2 , where (• • •) 
represents the temporal average, s is the location, and t 
is the time parameter. 

A link between each pair of sites si and S 2 on year y 
is defined based on the lagged cross-correlation function 
Xl, S2 (r > 0) = (e“(t)0? 2 (f + r)), and = 

X% (—r), where r is the time lag parameter that is 

considered up to 200 days (the results are not sensitive 
to this parameter). We also define the optimal time lag, 
t* , at which X% s (r) is maximal (or minimal), as the 
time delay of a pair Si, S 2 . When si is to the west of S 2 
and the time lag is positive, the link direction is to the 
east. We distinguish between positive and negative link 
weights as follows 


wv ,+ max(X| iS2 ) - meanpq^J 

Sl ’ S2 std(XlU) 


( 1 ) 


and, 


W v ~ 

VV Sl,S2 


min(Aj' i , S2 ) - mean (X^) 
std(X% uaa ) 


( 2 ) 


where max and min are the maximum and minimum 
values of the cross-correlation X% (r), mean and std 
represents the mean and standard deviation. Each cli¬ 
mate network is constructed based on 1 yr SSH data, and 
the time resolution for two consecutive climate network 
is 30 days. To eliminate the effect of slow modulations 
on correlation estimates, a two years high-pass filter has 
been used. Examples of SSH time series and their cross 
correlation-function is presented in Fig. SI of the SI. 

We construct climate networks for three sub-regions 
of the tropical Pacific: equatorial Pacific (2°S to 2°N), 
northern tropical Pacific (4°N to 7°N) and southern trop¬ 
ical Pacific (4°S to 7°S). To identify a threshold below 
which link weights W can be disregarded, we apply a 
shuffling procedure in which the order of years is shuf¬ 
fled while the order of days within each year remains 


unchanged. Thus, statistical properties of the data such 
as the distribution of values and their autocorrelation 
functions (within 1 year) are not affected by the shuffling 
procedure, but the statistical dependence between dif¬ 
ferent nodes is diminished. Fig. [lji depicts the temporal 
average positive link weights, (W + ), versus the geograph¬ 
ical distances of the links D , for the real and shuffled 
data in the equatorial Pacific regions (from 2°S to 2°N). 
High positive average link weight values which exist in 
the real data but not in the shuffled data are less likely 
to occur by chance. We thus consider only links that are 
separated from the shuffled links (indicated by the hor¬ 
izontal dashed lines and vertical arrows in Fig. [TJi); in 
addition, short distanced links (indicated by the vertical 
dashed line and the horizontal arrow in Fig. la) exhibit 
trivial correlations and are thus ignored (the “proximity 
effect” [H|). The probability density function (PDF) of 
average link weight in this region can be seen in Fig. S2 
of the SI. We find that significant links with a p-value 
p < 0.02 include 26.6% of the weighted links in the equa¬ 
torial Pacific (2°S to 2°N), 6.5% for 4°N to 7°N, and 
11.7% for 4°S to 7°S. 


One of the main characteristics of equatorial Kelvin 
waves is their phase speed (which is equal to their group 
velocity since these waves are non-dispersive)-below we 
use the links’ statistics to estimate this velocity. As¬ 
suming that the delay time r* is a good estimator for a 
dynamical interaction time between two nodes [22}, the 
phase velocities can be estimated by dividing the geo¬ 
graphical distance D between the nodes by the delay time 
r*, i.e., v ss D/t*. The PDF of phase speeds (of links 
with D > 500 km and p < 0.02) is shown in Fig. [jj. For 
the equatorial Pacific, the phase speed is positive (east¬ 
ward) and peaks around 3m/s, in agreement with the 
value and direction of Kelvin waves [§, [23|. 

The topological properties of the climate network in 
equatorial Pacific region show distinct behavior during 
El-Niho and normal times. The climate network consid¬ 
ered above is directed—positive and negative r* indicates 
eastward and westward flow respectively. We distinguish 
between a link that is pointing towards a node, or away 
from the node; the former/ latter is related to the total 
link weight pointing to/ from a specific node and is re¬ 
ferred to below as weighted “in”/ “out”-degree. Fig. |TJ>,d 
depict the “in” and “out” weighted degree versus longi¬ 
tude, averaging over all the latitudes from 2°S to 2°N 
during El-Niho and normal times, respectively. First, 
we find that the “out”-weighted degree is peaked around 
170°E at the equator with a much more pronounced 
peak during El-Niho, which indicates enhanced Kelvin 
waves activity, consistent with previous, more conven¬ 
tional, analysis [fj. Indeed, stronger Kelvin wave activity 
during El-Niho around 170°E was associated with strong 
WWBs [24| . which is enhanced wind activity that un¬ 
derlies the generation of equatorial Kelvin waves. Sec¬ 
ond, there is a peak in the “in”-weighted degree curve at 
~220°E during both El-Niho and normal times, where 
for El-Niho times the “in” degree curve is high all the 
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Velocity (m/s) Longitude (degree East) 

FIG. 1: Equatorial Pacific network (2°N to 2°S) analysis, (a) The temporal mean weight of positive links (W + ) versus their 
distance for both real (blue dots) and shuffled (red dots) data, (b) The “in” (red circles) and “out” (blue squares) weighted 
degree as a function of longitude for El-Nino times, for positive links (p < 0.02 and D > 500km). (c) The probability density 
function (PDF) of phase speeds that are based on positive links with p < 0.02 and distances larger than 500km. (d) Same as 
(b) but for normal times. 


way to the coast of South America while during normal 
times the “in”-degree curve decays for longitudes larger 
than ~220°E. We interpret this decay as strong dissipa¬ 
tion of the Kelvin waves, and our results suggest that 
the dissipation of Kelvin waves is much weaker during 
El-Nino. 

The attenuation of the Kelvin waves to the east of 
~220°E during El-Nino was implied based on wave pack¬ 
ets (found using Hovmollcr diagrams) using observation 
and model data [25], [26|. This attenuation may be at¬ 
tributed to several mechanisms, including the strong 
vertical current shear during La-Nina, strong wind at 
~220°E [25}, scattering of waves, and partial reflection 
of Kelvin waves into off-equatorial Rossby waves [2f|. 
Here, we use the climate network approach to quantify 
the waves properties in a statistical and consistent way. 

We next apply the climate network methodology to the 
northern and southern tropical Pacific to study Rossby 
(and other) tropical waves, respectively. The statistical 
properties of positive links in the climate networks in 
the northern tropical Pacific (4°N to 7°N) and southern 
tropical Pacific (4°S to 7°S), are shown in Fig. [2] The 
phase speed that is calculated based on links above signif¬ 
icance level [p < 0.02 and D > 500km) is negative (west¬ 
ward) and their PDF peaks around 0.5m/s and 0.6m/s 
respectively (Fig. [2L,d), in general agreement with the 
phase velocity of tropical instability waves (TIW) and 
off-equatorial Rossby waves 0 fi l27H30| . The structure 
of link weight versus distance is different for the north¬ 


ern and southern tropical Pacific (Fig. [5Ji,b); see [j] and 
supporting information (Fig. S3 and Fig. S4). Net¬ 
work analysis based on northern tropical Pacific is pre¬ 
sented in Fig. [3] The westward propagation of waves 
is evident during both normal and El-Nino times, with 
most of the wave activity concentrated between 200°E 
to 250°E (see also 0 )- Evidently, extra-tropical waves 
are less pronounced during El-Nino times (Fig. [3]). The 
enhanced wave activity during normal times can be as¬ 
sociated with the TIW. The characteristic phase speed, 
location and wave activity during El-Nino of TIW indeed 
fits the network observations [29|]. 

The TIW are not directly related to the theory of El- 
Nino even though they are possible sources of random 
perturbations that affect both the mean state and in¬ 
terannual climate variations. Therefore, we filter them 
out by employing a sixty days low-pass filter on the SSH 
data as their characteristic period is about 30 days, which 
is much smaller than the period of Rossby waves [30 i]. 
After this filtering, we reconstruct the climate network 
and calculate the corresponding wave characteristics (see 
Figs. S5, S6 of the SI). We relate this “filtered” cli¬ 
mate network with Rossby waves, which play an essen¬ 
tial role in the self-sustained El-Nino dynamics. Fig. Q] 
depicts the weighted degree at 4°N during El-Nino and 
normal times. The Rossby waves signal, represented here 
by the weighted degree of the climate network, is only 
clearly observed mainly to the west of ~ 240°E (120°W) 
(see Fig. [I]), although some of Rossby waves close to the 
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FIG. 2: The temporal average link weight (W + ) versus geo¬ 
graphical distance D in (a) the northern tropical Pacific (4°N 
to 7°N) and (b) the southern tropical Pacific (7°S to 4°S). The 
PDF of phase velocity in (c) the northern tropical Pacific and 
(d) the southern tropical Pacific. 


proach is useful to capture and quantify oceanic wave 
activity. This conclusion is further strengthened by the 
observation that the location of hubs in network coin¬ 
cide with the location of initiation and dissipation regions 
of equatorial oceanic waves. Moreover, the climate net¬ 
work shows distinct features during El-Niho and normal 
times, such as the pronounced equatorial Kelvin waves, 




FIG. 3: The out- (blue squares) and in- (red circles) weighted 
degree of nodes in the climate network from 4°N to 7°N during 
(a) El-Nino times and (b) normal times. 

eastern boundary (~ 280°E) are the reflection of Kelvin 
waves. Note that Rossby waves are more pronounced 
during El-Niho times: Around the date line (180°E) there 
is a peak in the weighted degree, which can be associated 
to the strong WWBs activity in this region. In addition, 
there is another peak in the weighted degree distribu¬ 
tion, between ~200°E to ~230°E, which can be associ¬ 
ated with the enhanced wind-stress curl in this region. 
This supports the recent suggestion that the enhanced 
wind-stress curl is related to the Rossby wave activity 0 ■ 
In summary, we construct and analyze the climate net¬ 
work of the tropical Pacific ocean based on daily satellite 
sea-surface height data (altimetry). We find that the 
most dominant long distance links in the network share 
the same characteristics (speed and direction) as equato¬ 
rial Kelvin waves, off-equatorial Rossby waves, and tropi¬ 
cal instability waves, indicating that climate network ap- 


FIG. 4: The weighted degree in 4°N for El-Nino (blue circles) 
and normal (red boxes) times. 

off-equatorial Rossby waves, and suppressed TIW dur¬ 
ing El-Niho. In addition, the dissipation of Kelvin waves 
are much weaker during El-Niho. This newly developed 
methodology based on long-term data (18 years) could 
serve as a quantitative tool to improve our understand¬ 
ing of the underlying mechanisms of climate phenomena 
like El-Niho. 

Recent studies have demonstrated the ability of 
climate network based-techniques to predict El-Niho 
events [H, [53] . This technique is based on atmospheric 
data above the sea surface, although as shown here the 
ocean plays a central role in El-Niho phenomenon. Based 
on the present study we conjecture that the inclusion of 
oceanic data into network analysis will improve the pre¬ 
diction power of network techniques. 
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